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Abstract. We explore the observational appearance of the photosphere of an ultrarelativistic flow with internal dissipation of 
energy as predicted by the magnetic reconnection model. Previous study of the radiative transfer in the photospheric r egion 
has shown that gradual dissipation of energy results in a hot photosphere. There, inverse Compton scattering of the thermal 
radiation advected with the flow leads to powerful photospheric emission with spectral properties close to those of the observed 
prompt GRB emission. Here, we build on that study by calculating the spectra for a large range of the characteristics of the 
flow. An accurate fitting formula is given that provides the photospheric spectral energy distribution in the ~ 10 keV to ~10 
MeV energy range (in the central engine frame) as a function of the basic physical parameters of the flow. It facilitates the 
direct comparison of the model predictions with observations, including the variability properties of the lightcurves. We verify 
that the model naturally accounts for the observed clustering in peak energies of the E ■ f(E) spectrum. In this model, the 
Amati relation indicates a tendency for the most luminous bursts to have more energy per baryon. If this tendency also holds 
for individual GRB pulses, the model predicts the observed narrowing of the width of pulses with increasing photon energy. 

Key words: Gamma rays: bursts - radiation mechanisms: general - methods: statistical 
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1. Introduction 

Despite the wealth of observational data obtained over the 
last few decades, the radiative mechanisms responsible for the 
prompt emission of a y-ray burst (hereafter GRB) remain elu- 
sive. From observations of thousands of bursts it has been 
shown that the GRB prompt spectral energy distribution is 
characterized by a peak (in E-f{E) representation) that clusters 
in the sub-MeV energy range. The peak smoothly connects to 
low and high frequency power-laws (Band et al. 1993; Preece 
et al. 1998; Kaneko et al. 2006). Theoretical arguments make 
use of the fact that a large fraction of the observed radiation is 
at photon energies above the pair creation threshold. Together 
with the huge luminosities involved in the GRB phenomenon, 
this implies that the emitting material is ultrarelativistic with 
bulk Lorentz factor F > 100 (for a review of these arguments 
see, for example, Piran 2005). Models for the prompt GRB 
phase thus have to account for both the acceleration of the flow 
and its observed spectral properties. 

For a flow to be accelerated to high bulk Lorentz factors, 
it must start with high energy-to-rest-mass ratio. Depending 
on whether the energy is in thermal or magnetic form, one 
has a fireball (Paczynski 1986; Goodman 1986) or a Poynting- 
flux dominated flow (Thompson 1994; Meszaros & Rees 1997; 
Spruit et al. 2001; Drenkhahn & Spruit 2002; Lyutikov & 
Blandford 2003; Uzdensky & MacFadyen 2006). In the fireball 
model, the flow goes through an initial phase of rapid acceler- 
ation where F ~ r until most of the internal energy has been 
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used to accelerate the flow (unless the baryon loading is very 
low and radiation decouples from matter before the accelera- 
tion phase is over). Further out, internal shocks that are a result 
of fast shells catching up with slower ones (Rees & Meszaros 
1994; Sari & Piran 1997) can dissipate part of the kinetic en- 
ergy and power the prompt emission. In magnetic models for 
GRBs, the flow acceleration is more gradual (e.g. Drenkhahn 
& Spruit 2002). 

Dissipation of magnetic energy through magnetic re- 
connection (Drenkhahn 2002) or current driven instabilities 
(Lyutikov & Blandford 2003; Giannios & Spruit 2006) can di- 
rectly power not only the prompt emission, but also the con- 
version of Poynting flux to kinetic energy of outflow. In fact, 
at high Lorentz factors magnetic dissipation is the only effi- 
cient mechanism for the conversion of Poynting flux to kinetic 
energy (Drenkhahn 2002, Drenkhahn & Spruit 2002, Spruit & 
Drenkhahn 2004; for a disagreeing view see Vlahakis & Konigl 
2003). This is discussed again in section 12". 1.1 1 

Unlike the internal shock model, where dissipation of flow 
energy is due to the variability of the central engine, variability 
and dissipation are, a priori, unrelated in magnetic dissipation 
models; dissipation takes place even in a steady magnetic out- 
flow. 

If the dissipated energy leads to fast particles, then the syn- 
chrotron or synchrotron-self-Compton mechanism appears as a 
natural one for the prompt emission. On the other hand, the 
large number of bursts with low energy slopes steeper than 
those that the synchrotron model can explain (e.g. Crider et 
al. 1997; Frontera et al. 2000; Girlanda et al. 2003), have trig- 
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gered alternative suggestions for the origin of the GRB emis- 
sion. These include saturated Comptonization (Liang 1997) 
and Comptonization by thermal electrons (GhiselUni & Celotti 
1999). 

Another possibility for the origin of the prompt emission 
is the photospheric emission of the flow (Meszaros & Rees 
2000; Ryde 2004; Thompson et al. 2006). A common feature of 
fireball and magnetic models is that they are characterized by 
(weaker or stronger) emission from the region where the flow 
becomes optically thin to Thomson scattering. If no energy dis- 
sipation takes place in that region, as in the case of fireball mod- 
els, the emission is expected to be quasi-thermal and a corre- 
sponding peak should be recognizable in the spectrum (Daigne 
& Mochkovitch 2002), contrary to observation. On the other 
hand, if there is substantial energy dissipation (through inter- 
nal shocks, magnetic reconnection or MHD instabihties) close 
to the photospheric region, the emitted spectrum can strongly 
deviate from a black body (Rees & Meszaros 2000; Pe'er et al. 
2005; Ramirez-Ruiz 2005). 

Giannios (2006, hereafter G06) studied the spectrum of 
the photospheric emission in the magnetic reconnection model 
(Drenkhahn 2002; Drenkhahn & Spruit 2002). The reconnec- 
tion model makes definite predictions for the characteristics of 
the flow and the rate of magnetic energy dissipation at diff'er- 
ent radii, enabling detailed study of the radiative transfer. This 
study showed that while radiation and matter are in approxi- 
mate thermodynamic equihbrium deep inside the flow, direct 
heating of the electrons through energy dissipation leads to 
an increase of their temperature even at optical depths ~ 50. 
Inverse Compton scattering of the underlying thermal radiation 
leads to a spectrum with a non-thermal appearance, peaking at 
~ 1 MeV (in the central engine frame) followed by a nearly flat 
high energy spectrum {E ■ f{E) ~ E^) extending to at least a 
few hundred MeV. 

Here, we build on that study with a wide parametric explo- 
ration of the spectral energy distribution of the photospheric 
emission in this model. We find that the photospheric spectra, 
computed from detailed Monte-Carlo simulations, can be ac- 
curately modeled by a smoothly broken power law model in 
the 10-10^ keV range. These are then used to derive analytic 
expressions for the dependence of the spectra on the character- 
istics of the flow. Since modulations of the characteristics of the 
flow by the central engine result in variable photospheric emis- 
sion, modeling of the spectral and temporal properties of the 
photospheric emission ofl'ers a unique probe into the operation 
of the central engine. 

We look at the properties of the distribution of the peak 
energy of the spectrum. The Amati relation is reproduced pro- 
vided that the more luminous flows are characterized by more 
energy per baryon (i.e. less baryon loading). If the same baryon 
loading-luminosity corelation holds during the temporal evolu- 
tion of GRBs, the model predicts the observed energy depen- 
dence of the width of the pulses of the GRB hghtcurves. 

In the next section, we review the main features of the mag- 
netic reconnection model and the treatment of the radiative 
transfer in the photospheric region. The analytical modeling of 
the photospheric spectra is given in Sect. 3 . Inferences concern- 



ing the properties of the central engine are given in Section 4 
and we summarize our findings in the Discussion Section. 

2. Dissipative GRB piiotosphiere 

The GRB flow, whether magnetized or not, has a characteris- 
tic radius at which it becomes transparent. This is the radius 
where radiation and particles decouple and the photospheric 
emission originates. If no energy is dissipated close to the pho- 
tosphere, the photospheric emission is expected to be quasi- 
thermal. On the other hand, significant dissipation of energy 
can strongly modify the photospheric spectrum. The physical 
source of energy dissipation can be the collision of shells of 
varying bulk Lorentz factors resulting in internal shocks in the 
flow (if, for some reason, the shells collide close to the pho- 
tosphere). Dissipation near the photosphere is generic in the 
magnetic reconnection model, where the dissipation of mag- 
netic energy is gradual and takes place over several decades of 
radius, typically including the photosphere. 

Unlike the internal shock model, where dissipation of flow 
energy is due to the variability of the central engine, variability 
and dissipation are, a priori, unrelated in the magnetic recon- 
nection model; dissipation takes place even in a steady mag- 
netic outflow. For typical inferred GRB parameters, the dissi- 
pation in a magnetic reconnection model takes place around the 
photosphere of the flow (Drenkhahn 2002). The time for the 
flow to reach the photosphere as observed by an on-axis exter- 
nal observer is so short that the flow can be treated as approxi- 
mately steady for all but the shortest observed time scales. This 
greatly simplifies the interpretation of light curves and varia- 
tions in the energy spectra. Of course, this applies only if it can 
be demonstrated that a substantial fraction of the emission in 
the 10-1000 keV range is in fact produced around the photo- 
sphere of the flow. The spectra obtained in G06 are promising 
in this respect. 

The reconnection model makes specific predictions for the 
rate of energy dissipated as a function of radius; this makes 
it possible to study the radiative transfer in the flow in detail 
(G06). Analytic estimates demonstrate that the flow develops 
a hot photosphere where inverse Comptonization results in up- 
scattering of the photons advected with the flow. The numerical 
(Monte Carlo) calculations in G06 show that the spectra have 
many properties of the observed prompt GRB emission, for a 
wide range of baryon loading and luminosities of the flow. 

In the following subsections, we review the basic charac- 
teristics of the recormection model and of the radiative transfer 
in the flow. 

2.1. The reconnection modei 

An important physical quantity in the reconnection model is the 
ratio ctq of the Poynting flux to kinetic energy flux at the start of 
the flow, the Alfven point vq. This is a quantity equivalent to the 
'baryon loading' parameter in fireball models. It determines the 
terminal bulk Lorentz factor of the flow Too ~ <^q^^- The flow 
must start Poynting-flux dominated with ctq > 30 for it to be 
accelerated to ultrarelativistic speeds with Too > 100 that are 
relevant for GRB flows. 
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In the reconnection model, the magnetic field in the flow 
changes polarity on a scale A. If the magnetic field anchored 
in the rotating central engine is nonaxisymmetric (the 'AC 
model'), this scale is (in the central engine frame) of the order 
of the light cylinder rj, (A ^ ro ^ 2jtc/Q. = 2nr]) where O is the 
angular frequency of the rotator This is as in the oblique rotator 
model for pulsar winds (Coroniti 1990). The rate of magnetic 
reconnection in the model is parameterized through the veloc- 
ity Vrec with which magnetic fields of opposite direction merge. 
As in most models of magnetic reconnection, v,ec scales with 
the Alfven speed, va, Vrec = sva- A nominal value used for e 
is 0.1. For the flows with cro » 1 that are of interest here, the 
energy density of the magnetic field is larger than the rest mass 
energy density, hence va ~ c, and the reconnection takes place 
with subrelativistic speeds. For a detailed analysis of relativis- 
tic reconnection see Lyubarsky (2005). 

2.1 .1 . Poynting flux conversion 

If the energy extracted from the central engine is initially in 
the form of a Poynting flux, acceleration of the flow requires 
efficient conversion of this flux into kinetic energy. The stan- 
dard mechanism for acceleration of nonrelativistic flows, mag- 
netocentrifugal acceleration, (Blandford & Payne 1982) allows 
conversion of a substantial fraction (some 40%), with the rest 
remaining in the flow as magnetic energy. For relativistic flows, 
the conversion efficiency for this same mechanism is much 
lower (Michel 1969), so that in the case of a GRB almost all 
energy would end up in the afterglow rather than as prompt 
emission. 

Arbitrarily high effiency of conversion is possible, however, 
if there is a channel by which the magnetic energy carried by 
the flow can dissipate by internal processes. Such dissipation 
results is an outward decrease of magnetic pressure which ac- 
celerates the flow efficiently (Drenkhahn 2002; for a more de- 
tailed physical explanation of this at first sight somewhat coun- 
terintuitive process see Spruit & Drenkhahn 2004). 

The dissipation is gradual and takes place up to the 'sat- 
uration radius' r^, the distance where reconnection stops and 
the flow achieves its terminal Lorentz factor. During the mag- 
netic dissipation process, about half of the Poynting flux di- 
rectly converts into kinetic flux (leading to acceleration) and 
the other half into internal (thermal) energy of the flow. 

If the dissipation takes place under optically thin condi- 
tions, the thermal energy can be radiated away quickly. In this 
case half of the Poynting flux is converted to radiation, the other 
half to kinetic energy of the bulk flow. If radiation is inefficient, 
on the other hand, the internal energy adds to kinetic energy of 
the bulk flow, through adiabatic expansion (as in the case of a 
pure hydrodynamic fireball). 

2.1 .2. Properties of the magnetically accelerated flow 

The reconnection model predicts gradual a acceleration of the 
flow r ~ r^^^ in the regime ro r r^, while no further accel- 
eration takes place beyond the saturation radius. Summarizing 
the results derived in Drenkhahn (2002), the bulk Lorentz fac- 



tor of the flow is approximately given by 
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In the steady spherical flow under consideration the comov 
ing number density can be written as 

L 



(3) 



where L is the luminosity per steradian of the GRB model. 
(In this form the expression can be compared with the fireball 
model, where the baryon loading parameter t] replaces the fac- 
tor crl^^). With eqs. O, (O we have 
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The comoving magnetic field strength B' is given by 



47rL 1.4-10** L, 
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G, for r < r^ (5) 



with the magnetic field in the central engine frame given by 
B - TB'. [The usual notation A - WA^ is used; the 'reference 
values' of the model parameters are cro = 100, s - 0.1, Q - 
10'^ rad-s-', L = lO^^ erg-s-'-sterad"'.] 

2.2. Radiative transfer in the photospheric region 

In addition to the saturation radius rs, another characteristic ra- 
dius of the flow is the Thomson photosphere. A photon moving 
radially outward in a flow with F » 1 encounters an optical 
depth dr = n'crTdr/2F (Abramowicz et al. 1991). Integrating 
from r to 00 the characteristic Thomson optical depth as a func- 
tion of radius is 

20 L52 
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Setting T = 1 yields the photospheric radius r. 
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In deriving these expressions, we have assumed that rph < rj. 
At high baryon loading, the Lorentz factors can become so low 
that most of the Poynting flux dissipates akeady inside the pho- 
tosphere, i.e. rph > r,. With the expression for the number den- 
sity valid beyond the saturation radius from (|4|, a similar cal- 
culation gives the radius of the photosphere in this case. 
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Due to the dominance of scattering, the details of radia- 
tive transfer become important already at fairly large optical 
depth in the flow. Equilibrium between radiation and matter 
holds only at Thomson depths greater than about 50. At optical 
depths 0.1 < T < 50 the electron distribution stays thermalized, 
but is out of equilibrium with the photon field. Compton scat- 
tering of the photons needs to be treated in detail in this region. 
This has been done in G06, using both analytical and numerical 
tools. One of the main findings of this study is that the nature 
of the photospheric spectrum is determined by the location of 
the photosphere relative to the saturation radius. If rph » r,, all 
the energy dissipation takes place in optically thick conditions, 
the radiation is efficiently thermalized and suffers substantial 
adiabatic losses. On the other hand, if r, > rph, energy dissipa- 
tion at moderate and low optical depths leads to a photospheric 
emission that has a highly non-thermal appearance. 

There is thus a characteristic spectral transition of the pho- 
tospheric emission depending on the relative locations of the 
photospheric and saturation radii. In terms of the physical prop- 
erties of the flow, it takes places at the critical value of the 
baryon loading cro,cr, for which rph - r^. For cr > cro,cr, '"ph < '"s- 
Using eqs. d?) and ^ one finds 
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At radii much shorter than the photospheric radius, the re- 
leased energy is efficiently thermalized and shared between 
particles and radiation. In this region, the comoving tempera- 
ture Tth of the flow is calculated, under the assumption of com- 
plete thermalization, by integrating the energy released at dif- 
ferent radii in the flow and taking into account adiabatic cool- 
ing. This leads to (Giannios & Spruit 2005) 
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Complete thermalization of the dissipated energy 
achieved up to a radius defined as the equilibrium radius 
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where /e stands for the fraction of the dissipated energy that 
heats the electrons (rather than the ions) and is assumed to 
be of order of unity. As a reference value we adopt /e - 0.5. 
Note that the equilibrium radius lies at a distance one order of 
magnitude below the Thomson photosphere. At radii r > rgq, 
direct electron heating through magnetic dissipation leads to 
electrons with characteristic energy higher than that of the ra- 
diation field and upscattering of the photons takes place in this 
region. Because of the rather steep increase of the electron tem- 
perature with radius above rgq, this upscattering becomes effi- 
cient especially around the location of the photosphere, where 
the comoving electron temperature is of the order a few tens 
keV. For these temperatures and optical depths the result is a 
hard power law tail in the spectrum, characteristic of unsatu- 
rated Comptonization. 

For the detailed calculation of the emergent spectrum and 
the electron temperature above the equilibrium radius we use 
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Fig. 1. Photospheric spectrum of the reconnection model for 
GRB flows (photon energies in the central engine frame). For 
moderate and high initial Poynting flux ratio cro the spectra 
peak around 1 MeV. The flat high energy tail is the result of 
inverse Compton scattering of photons off hot electrons close 
to the photospheric radius. For low o-q dissipation stops below 
the photosphere, and the resulting spectrum is quasi-thermal. 



a Monte Carlo method as described in G06. It treats Compton 
scattering in a medium with density and bulk Lorentz factor 
given by eqs. (|4]l and ([U respectively. The special relativistic 
effects related to both the bulk motion and the scattering cross- 
section of photons in the flow are taken into account. The inner 
boundary of the computational domain is the equilibrium ra- 
dius defined by eq. (fTOl i. There, a black body photon spectrum 
is injected with temperature given by eq. (|9]l. The outer bound- 
ary rout is, rather more arbitrarily taken at a radius where the 
optical depth becomes small (the value t = 0.1 is used for the 
rest of the Section, but see Section ( 13.1b where the dependence 
of the emitted spectrum on the choice of rout is investigated). 

The electron temperature as a function of distance in the 
flow is determined from the balance of radiative cooling by 
Comptonization and heating of the electrons by magnetic dissi- 
pation. The numerical results verify the analytical expectation 
that the electron temperature increases rapidly in the photo- 
spheric region and that significant Comptonization takes place 
there (G06). 

A selection of characteristic emergent spectra from the cal- 
culations are presented in Fig.[T]in E ■ f{E) representation for 
different values of the parameters of the model. The spectra are 
plotted in the central engine frame with arbitrary normalization 
(the energetics are discussed in the next Section). The peak of 
the spectrum is close to 1 MeV for a large region of the param- 
eter space which is relevant for GRB flows. The dependence of 
the peak energy on the parameters of the model is very weak. 
This can be understood from the fact that the peak energy is re- 
lated to the peak of the thermal emission advected with the flow 
at the equilibrium radius which, in the central engine frame, is 
£'p,th <^ r(req)rth(req) oc Lj^^(eQ)3^ Vq^'j , where in the 

last step the expressions ©, (HJ and (fTOl i have been used. The 
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Fig. 2. Smoothly-broken power law fits (dashed) to the spectra 
of Fig. [T] (solid). Photon energies in the central engine frame. 
The 5-parameter SBPL formula (eq. (fTTl i) fits the numerically 
calculated spectra very accurately in the 10 keV to 10 MeV 
energy range. 

detailed numerical investigation leads to very similar depen- 
dence of on the model parameters (see eq.fTsll. 

The inverse Compton scattering that takes place above the 
equilibrium radius leads to a moderate increase of the energy 
of the peak and, most important, to a power-law high energy 

emission with photon number index 2.5 that extends up to 

a few hundred MeV. This high energy part of the spectrum is 
the result of unsaturated Comptonization taking in the t < 1 
region and appears for flows with ctq ^ cro,ci . 

In the lower ctq cases, the emergent spectrum has much 
weaker emission above its peak. In these cases the magnetic 
dissipation stops below the Thomson photosphere and there 
is only weak Compton upscattering taking place in the photo- 
spheric region. Such quasi-thermal emission has been observed 
in a fraction of GRBs (Ryde 2004; 2005) suggesting these may 
be bursts of low ctq (high baryon loading). 

3. Analytical fits and parametric study of the 
spectra 

The spectra obtained with our Monte-Carlo radiative transfer 
calculations have systematic dependencies on the parameters 
of the model (initial Poynting flux to kinetic flux ratio cry, lu- 
minosity per steradian L and reconnection rate which appears 
in the combination ef2). In this section we condense these de- 
pendences into a more practically useful form by developing 
fitting formulas. 

After various attempts to find analytical functions that 
fit the numerically computed spectra, it turned out that the 
smoothly connected power law (hereafter SBPL model Preece 
et al. 1994; Ryde 1999; Kaneko et al. 2006) provides a very 
good description in the 10 keV to 10 MeV energy range, which 
is the most relevant for observations of the prompt emission 
of GRBs. The SBPL is a generalization of the sharply broken 
power law model by the addition of a parameter governing the 



width of the transition from the low energy spectral slope to 
the high energy slope. The fitting model has five parameters: an 
amplitude A, low and high photon indexes X\ and X^, a break 
energy /Sb and a break scale A. The transition in slope is de- 
scribed by a hyperbolic tangent function. 

The SBPL model fits the data (in our case the spectrum) 
with the expressiorQ (for the derivation see for example Ryde 
1999) 



I E e'' +e-i \- 
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ln(£/£b) 



A 



ln(100/£b) 



(11) 



(12) 



The amplitude A is in photons s ' and all energies in keV. 

The break energy Eh is related to the peak energy Ep of the 
distribution Ef{E) by 



A2+2 



A/2 



(13) 



(e.g., Kaneko et al. 200611 

Since this fitting model does not have a high energy expo- 
nential break, it cannot describe the spectra above -100 MeV. 
The SBPL fits also show small systematic deviations for ener- 
gies below -several keV. Thus, we limit our fits to the, never- 
theless broad, 10 KeV - 10 MeV energy range. We have used 
the SBPL model to fit a large number of theoretical spectra for 
models of varying luminosity L, baryon loading ctq and recon- 
nection rate parameterized by eQ. 

A few, rather representative, examples of such fits are 
shown in Fig. |2] The fits (shown as dashed curves) lie very 
close to the numerical results with the exception of the lowest 
end of the energy range considered. Overall, SBPL fits are a 
very good way to summarize the photospheric spectrum pre- 
dicted by the dissipation model. 

With these fits, the low and high energy slopes Ai and A2 
and the peak energy E^ of the spectrum can be presented as 
functions of the parameters of L, ctq and eQ of the model. 

Fig. |3]shows the ratio of the photospheric emission to total 
flow luminosity, i.e. the (photospheric) radiative efficiency of 
the model and Fig. |4]the peak energy Ep of the spectrum. For 
(To around the critical value o-Q a, the radiative efficiency peaks 
at w 0.35/e, independent of the luminosity of the flow L and the 
reconnection rate parameter sQ. . For o-q < ctq cr, the magnetic 
dissipation stops below the photosphere and radiation suffers 
substantial adiabatic cooling before it decouples from matter; 
this results in a steep increase of the radiative efficiency and of 
the peak energy as functions of cro. Changes in L and sQ. lead 
to characteristic displacements of the curves in the Lph/L - ctq 
and Ep - ctq planes (see Figs. 3 and 4). 



' We choose a representation of the SBPL model where only the 
natural logarithm appears. This is achieved by suitably renormalizing 
the break scale with respect to that used by Ryde (1999). 

- The very existence of a peak in the spectrum demands that /ii +2 > 
and A2 + 2 < 0, which almost always holds here. 
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Fig. 3. The (photospheric) radiative efficiency of the flow (ratio 
of photospheric luminosity to flow luminosity), as a function 
of the initial Poynting flux to kinetic energy flux ratio ctq. At 
low (To, most of the dissipation happens at great optical depth, 
and the photospheric radiation field is reduced by adiabatic ex- 
pansion. At high ctq most dissipation takes places beyond the 
photosphere. The solid red line shows the effect of increasing 
the electron fraction of the dissipated energy from 50 to 100%. 
Changes in luminosity L and eQ are equivalent to shifts along 
the cr-axis. The curves are fits to smoothly broken power laws 
(best fit values in the appendix A.2). 




Fig. 4. Dependence of peak energy (peak of the photospheric 
E ■ f{E) spectrum) on baryon loading, flow luminosity and re- 
connection speed. Above the characteristic Poynting flux pa- 
rameter ctq^, (Eq. 7), depends only weakly on baryon load- 
ing. At high baryon loading (low ctq), the spectrum becomes 
quasithermal (cf. Fig. 1), peaking at a low energy. Curves show 
the smoothly broken power law fits (cf. appendix A.2). 



Analytic fits for the radiative efficiency and peak energy as 
functions of the characteristics of the flow are given in the ap- 
pendix. These expressions reduce to simple power laws when 
~ cTcr. Thus, for a sufficiantly high initial magnetization pa- 



rameter in the flow, the radiative efficiency has the following 
dependence on the model parameters (see Sect. A.2) 

Lph/L ^ 0.07(roi/(eQ)°-2L02\ cr,, > cr„,er. (14) 

Similarly, for ctq > ctq ci the expression for the peak energy 
simplifies to (see Sect. A. 3) 

£p ~ 12004^2 '(eil)3 ''VJJj keV, o-q > cro.ci-. (15) 

The high-energy photon index A2 is almost constant at 
2.5 for (To > CTQ^cr and decreases rapidly for smaller ctq (be- 
cause there is little or no dissipation close to the photosphere 
for low values of ctq). The expression for A2 as a function of 
the characteristics of the flow is given in the appendix A.4. The 
low-energy photon index A\ shows very small variations and 
remains approximately ~ in most the parameter space inves- 
tigated. For simplicity, we keep Ai - Q independently of the 
characteristics of the flow. On the other hand, significant soft- 
ening of A\ appears when the spectral peak approaches the low- 
est energy range over which the spectral fits are performed (i.e. 
at 10 keV). This may lead to deviations of the fitting formu- 
las and the actual predictions of the photospheric model in the 
~ 10 keV range (in the frame of the central engine). The break 
scale A is practically constant -1.6 for high ctq and makes a 
transition to the value ~ 2.8 for cry < cro cr. This transition can 
be modeled by a hyperbolic tangent (see appendix A.2). 

3.1. The non-photospheric emission 

The spectra above were obtained with an outer boundary set at 
a scattering optical depth of 0. 1 . In the optically thin part of the 
flow (t < 0.1) the electron temperature will continue to rise if 
significant dissipation takes place there, and this will contribute 
to the Comptonization. We have not included this contribution 
so far because the assumption of a thermalized electron distri- 
bution becomes questionable at t<0. 1 , and a nonthermal distri- 
bution will change the Comptonization process. This makes the 
emission in these regions dependent on poorly known details of 
the magnetic dissipation process. 

Notwithstanding this, it is possible to get an idea concern- 
ing at which energy band the optically thin emission may be 
expected to take place by extending the outer boundary of 
the Monte Carlo calculation to larger radii (or equivalently to 
smaller optical depths in the flow). In this excercise we assume 
that the electrons still follow a thermal distribution in this re- 
gion. In Fig. (|5]l the emitted spectrum for the reference values 
of the parameters of the model is shown for various values of 
the optical depth Tout at the outer boundary of the numerical 
calculation. The spectrum below the peak energy remains es- 
sentially the same. The higher energy emission is significantly 
enhanced as a result of scattering in the very hot, optically thin 
part of the flow. The high-energy spectral slope hardens, ap- 
proaching A2 ~ -2, and there is substantial emission in the 
~GeV energy range. 

The vertical dashed lines show the energy range used for 
the fitting formulas presented in the previous section. The op- 
tically thin emission affects the high-energy spectral slope A2 
somewhat, which now lies in the 2 to -2.5 range. The opti- 
cally thin emission increases the luminosity in the 10-lOOOkeV 
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Fig. 5. Modification of the spectrum when an estimate of the 
emission in the optically thin part of the flow is included. 
[Electron distribution is assumed to remain thermal at low opti- 
cal depths.] Parameter is the optical depth at the outer boundary 
of the calculation (photon energies in the central engine frame). 
Dashed lines mark the range used for the fitting formulas. 

range only moderately, but can produce a significant peak in 
the GeV range. 

4. Interpretation of the observations 

In the following we use a first comparison with observations to 
make inferences about the properties of GRB outflows in the 
magnetic dissipation model. 

4.1. The value of Ep 

The observed spectral peak of the prompt emission spectrum of 
GRBs has a narrow distribution peaking at -250 keV (Band et 
al. 1993; Kaneko et al. 2006); by far the largest amount of data 
emanates from BATSE. To compare this distribution with the 
prediction of the photospheric model ( fTSl l the redshift distribu- 
tion of GRBs has to be taken into account. This distribution is 
not well known for the BATSE bursts, but various estimates in- 
dicate that an average value of z ~ 1 .5 - 2 is reasonable. For a 
typical burst at this redshift the model predicts that (cf. eq. (TSi ) 



pobs _ 



1 -HZ 



- 4504'i(eQ)'3"^ 



.0.3 

0,2 



keV. 



(16) 



The dependence on the characteristics of the flow is weak; 
the model thus predicts a narrow distribution of the peak en- 
ergy of bursts, in agreement with observation. For the reference 
values of the parameters used in eq. (fTST l. the peak energy pre- 
dicted by the model is at 450 keV, a factor of ~ 2 above the 
value for BATSE bursts. Though this is not a large factor, the 
weak dependence on model parameters implies a significant 
difference in one or more of them compared with the reference 
value. A low value for ctq <k 100 is not attractive, since then the 
flow becomes too slow to explain the typical GRB. The very 
weak dependence of the peak energy on the luminosity also 



excludes the burst luminosity. Since the central engine is most 
likely a stellar mass compact object, the rotation rate, Q, must 
be of the order 10^-10** rad s"' . The most uncertain parameter 
of the magnetic dissipation model is actually the reconnection 
speed parameter s, for which we have so far assumed a value 
~ 0.1 (Drenkhahn 2002). Q. = 10'* rad s ' is perhaps on the 
high side for the torus in a collapsar, and e = 0. 1 possibly too 
optimistic. A value (eQ) ~ 100 rad s~' instead of 10^, which 
would bring the predicted value of Ep into agreement with ob- 
servation, can easily be accommodated within these uncertain- 
ties. 

This indicates that magnetic dissipation timescale is one or- 
der of magnitude longer than that assumed by the reference 
values of the parameters. In the following we adopt (eQ) = 
10~(eQ)2 as a reference value for the parameterization of the 
dissipation timescale. 

4.2. The Amati relation 

There are indications in the literature that the peak energy Ep 
of a burst correlates with the isotropic equivalent energy: Ep oc 
E^^ (Amati et al. 2002; Amati 2006), although this correlation 
is probably not universal (Nakar& Piran 2005; Band & Preece 
2005). The isotropic equivalent energy of the prompt emission 
is Eiso ~ LphtcRB, where tcRs is the duration of the burst. For 
an estimate of the energy of the prompt emission, the relevant 
luminosity is that of the photosphere Lph and not that of the 
flow L. 

Using the expressions ( fTsT i and ( fT4b one can find the depen- 
dence on the parameters of the flow of the ratio 



e"' 

ISO 



GRB 



(eD.)' 



,0.23^1.1 



(17) 



0-, 







According to the Amati relation, this ratio stays approximately 

constant over a wide range of Eiso- Since the dispersion in the 

1 y 

distribution of durations of f^'^g for long GRBs is rather small 
~ 10'^^, it is unlikely that a dependence in a systematic way of 
fcRB with another parameter (say the luminosity) is responsi- 
ble for keeping the ratio ( fTTl ) constant. On the other hand, the 
ratio is approximately constant for ctq oc /,~0-4-o.5 different 
burst^. The Amati relation, therefore, indicates that there is a 
corelation of baryon loading and luminosity, in the sense that 
the most luminous GRB flows have more energy per baryon. 
We see in the next Section that if the same baryon loading- 
luminosity corelation holds during the evolution of individual 
GRBs, the model predicts the observed energy dependence of 
the width of pulses in the GRB lightcurves. 

4.3. GRB variability and spectral evolution 

Up this point, we have dealt with a steady flow that is charac- 
terized by fixed luminosity and baryon loading. This approach 
can, to some extent, describe the rough, time average, prop- 
erties of a GRB. On the other hand, the GRB lightcurves are 



^ Assuming that (eSI) does not change much from burst to burst; if 
it varies in a systematic way with, say, the luminosity, the scaling of 
(Tq with L can be somewhat different. 
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highly variable and are composed by a sequence of pulses be- 
lieved to be the "building blocks" of their lightcurves (e.g. 
Fenimore et al. 1995). The question that rises as to where 
the variability comes from in the photospheric model for the 
prompt emission and under what conditions a steady state 
model such as the magnetic reconnection one can describe a 
variable flow. 

An ultrarelativistic flow that varies on a timescale 6t can be 
still described by a steady state model as long as 5t > r/F^c. 
When the last condition holds, the leading and trailing parts of 
a piece of the flow with width 6r ^ c6t have not come into 
causal contact. As long as this is the case, different parts of 
the flow behave as part of a steady wind. A robust proof of 
this statement can be found in Piran et al. (1993) for fireballs 
and Vlahakis & Konigl (2003) for MHD flows. The relevant 
radius of interest in this study is the photospheric radius. As 
far as the photospheric emission is concerned, the steady state 
approximation is valid for a flow which varies on timescales 



6t > 



rph 



r2(rph)c 



2 ■ 10" 



rl/5 
-"52 



s, 



(18) 



0,2 



where eqs. ^ and ([T]l are used in the last step. The flow can 
treated as approximately steady for all but the shortest observed 
(sub-millisecond) timescales. 

4.3.1 . Central engine: the stronger the cleaner 

Variability and dissipation are, a priori, unrelated in the mag- 
netic reconnection model in which dissipation takes place, even 
in a steady magnetic outflow. On the other hand, the flow does 
evidently evolve during a GRB. In the context of the present 
model, the observed variability reflects changes in the lumi- 
nosity (and possibly the initial magnetization parameter or even 
the reconnection rate eQ) during the burst. Since the flow up to 
the photosphere can be treated as quasi-stationary for all but the 
shortest time scales observed in a burst, the variation of spectral 
properties during a burst directly reflects variations in the cen- 
tral engine. This is in contrast to models in which the prompt 
radiation is produced at much larger distances from the source, 
such as external shock models. It is also in contrast with the 
internal shock model, since the internal evolution of the flow 
between the source and the level where radiation is emitted is 
a key ingredient in this model. Deducing properties of the cen- 
tral engine from observed burst properties is thus a much more 
direct prospect in the magnetic dissipation model. 

As a first application we consider here the width of sub- 
pulses in a burst. The observations show this width to be a de- 
creasing function of photon energy (e.g. Fenimore et al. 1995). 
GRBs with complex light curves (composed of many sub- 
pulses) follow the same Amati relation as simple ones (with 
one or a few pulses). In the above, we have shown how this 
relation implies a dependence of the baryon loading of the 
burst on its luminosity. We now make the logical extrapola- 
tion of assuming that this same relation also holds for the time- 
dependent flow during a burst. Varying the luminosity of the 
flow with time in such a way that it mimics the typical pulse 
profile, one can simulate the pulse appearance predicted by the 



bandl (25-55 keV) 
baii<i2 (55-115 keV) 
baiid3 (1 15-320 keV) 
bamU (320-2000 keV) 




Fig. 6. Narrowing of pulse width with photon energy, as pre- 
dicted by the model for individual subpulses of the burst. The 
assumption is that the relation between luminosity and baryon 
loading deduced from the Amati relation, <to oc L" "*, holds also 
within a burst. 

photospheric model in different energy bands. To reproduce the 
characteristic pulse shape, we consider time dependence of the 
luminosity of the flow of the form (cf. Norris et al. 1996) 



L(t) 



-\t/t,\ 



t < 



t > 



(19) 



where fr - 0.2/(1 -H z) s and fd = 0.5/(1 -H z) s are the rise 
and decay timescales of the lightcurve (in the central engine 
frame) respectively and L^ax = 10^^ erg s"' sterad"'. We also 
assume that the baryon loading tracks the luminosity with the 
expression ctq 2 = (indicated by the Amati relation) and set 
sQ. = 100 rad s The redshift of the burst, needed to correct 
for photon reshifting and cosmological time dilation, is taken 
to be z = 2. 

Having specified the parameteres of the flow as functions of 
time, one can calculate the lightcurve (in the observer's frame) 
of the photopsheric emission at different energy bands pre- 
dicted by the photospheric model using eq. ( fTTT i and the expres- 
sions in the appendix. In fig.|6l the predicted lightcurves in the 
BATSE energy bands are plotted. The lightcurves of the harder 
energy bands are narrower than those of the softer The width 
in the 320-2000 keV band being approximately half of that of 
the 25-55 keV band close to what is observed. The energy de- 
pendence of the lightcurves comes from the spectral hardening 
of the photospheric emission with increasing luminosity L (and 
therefore increasing photospheric luminosity Lph; see eqs. (fT4l i 
and (fTSIl). This spectral hardening results in relatively enhanced 
emission in the hard bands (with respect to the softer ones) 
when the luminosity is close to its peak. 

5. Discussion/conclusions 

In the magnetic reconnection model, the flow that powers a 
GRB is initially Poynting flux dominated. The model is char- 
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acterized by gradual dissipation of magnetic energy that leads 
to acceleration of the flow and an increase in its internal en- 
ergy. For characteristics of the flow typical to those expected 
for GRB flows (i.e., high luminosities and low baryon load- 
ings) the dissipation continues at distances outside the point 
where the flow becomes transparent with respect to Thomson 
scattering (i.e. outside the Thomson photosphere). The energy 
release close to the photospheric radius leads to a hot photo- 
sphere. Inverse Compton scattering of the advected with the 
flow photon field by hot electrons leads to photospheric spectra 
with a characteristic non-thermal appearance and with peak en- 
ergy and spectral slopes close to those observed in the prompt 
phase of GRBs. The radiative transfer problem has been stud- 
ied in detail, with both analytical and numerical tools, in G06. 

Here, we build further on that study by providing an accu- 
rate analytical description of the spectra predicted by the model 
in the x-ray and y-ray regime. We show that the photospheric 
emission from the reconnection model is well described by a 
smoothly broken power law model. Simple fitting formulas are 
presented for the bolometric luminosity and peak energy of the 
E ■ f{E) spectrum as a function of the flow parameters (i.e. the 
luminosity of the flow L, the baryon loading ctq and the recon- 
nection timescale parameter eQ). 

The model predicts the observed clustering of the peak en- 
ergy of the E ■ f{E) spectrum in the sub-MeV range and its 
quantitative comparison with observations indicates that the 
combination of parameters (eD) (where e stands for the frac- 
tion of the Alfven speed with which reconnection takes place 
and Q for the angular frequency of the rotator) has values 
(sD.) ~ 100 rad s ' . This is compatible with reconnection speed 
proceeding with e ~ 0.03 and a rotator with angular frequency 
Q ~ 10^-^rad s"'. The Amati relation (i.e. Er, ~ Ef'^ Amati 

^ ISO ' 

et al. 2002; 2006) can be explained by the model provided that 
the most luminous bursts are characterized by higher energy- 
to-rest-mass ratio. 

More constraints to both the model and the operation of the 
central engine can come from the study of the variability prop- 
erties of the GRBs. In the photospheric model, the observed 
variability is a direct manifestation of the activity of the central 
engine. Individual pulses, considered to be the building blocks 
of the GRB lightcurves (e.g. Fenimore et al. 1995) come from 
modulations of the luminosity and the baryon loading of the 
flow and are ultimately linked to the central engine variabil- 
ity. Modeling of the observed properties of the GRB pulses 
can, thus, provide us with a unique probe into oparation of the 
central engine on diferent timescales during the GRB prompt 
phase. 

We have simulated GRB pulses by assuming that the lumi- 
nosity of the flow varies following the typical observed pulse 
profile (cf. Norris et al. 1996) and that the magnetization pa- 
rameter (To corelates with the luminosity of the flow as indi- 
cated by our interpretation of the Amati relation. The pulse 
width predicted by the model is narrower in the harder energy 
bands, in agreement with observation. 
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Appendix A: Fitting functions for the spectral 
parameters 

Here we give the expressions that describe how the photo- 
spheric luminosity Lph, the peak energy E^, the high frequency 
power law and the spectral break A depend on the physical 
parameters L, cry and of the reconnection model. We first, 
derive some useful expressions related to the SBPL model that 
will simplify the apearance of the fitting formulas. 
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A.1. A note on the Smoothly connected power-law 
model 
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The fit for the peak energy of the E ■ f{E) spectrum on ctq 



IS 



The expression that describes the SBPL model can be quite ^ ' (e* + e"*) 

generally written (e.g., Ryde 1999) ""UOO/ I 100 / I / 



where 

ln(x/xb) 



A 



^piv 



ln(xpiv/xb) 
A ■ 



(A.l) 



(A.2) 



with 



Ap ^ 1200(eQ)5-^^L°2" 



keV. 



For (To > ctq^j this reduces to (cf. eq. ( |A.4| l) 



1200(sfi)°-33L«2"cr«;^ 



keV. 



The value Xpiv is a reference value of the x variable one has the 
freedom to choose. For x - jCpiv one easily check that gix^n) - 
A. 

This expression can be simplified if the reference value 
Xpiv is taken sufficiently far from the break x\,. If, for example, 
g9piv 3> g-?piv^ one can simplify the expression for the SBPL 
model. 

Under the condition that e'^p" » g-9piv eq. (lA.U can be 
brought into the form 

g{x)^Ai—] (^] ' ie'^+e-"] . (A.3) 



The fit for the high-energy slope is 



(A. 10) 



(A. 11) 



(A. 12) 



(A.13) 



The break scale is approximately constant, A ^ 1.6, for 
high (To and makes a transition to A ^ 2.8 for low values of (Tq. 
The transition is modeled with the hyperbolic tangent function 
with best fit 



A : 



2.2-0.6tanh(8^^i^^) 



(A. 14) 



An approximation to the spectrum itself, as a function of 
baryon loading, flow luminosity and reconnection rate, is given 
One can also derive the asymptotic expressions for the low and by eq. (HB, with parameters given by the fitting formulas 



high power-law dependence of g{x) on x in the limits e'? » e ■? 
and vice versa. So expression dA.SI l gives 



g(x)^Al—f, e^» 
and 



e ' «c e 



(A.4) 



(A.5) 



A^piV / 



The form (IA.3l l of the SBPL model is used below for fitting the 
dependency of the photospheric spectrum on the parameters of 
our magnetically powered GRB flow. 



above. The normalization factor A in eq. ( fTTI ) follows by equat- 
ing the integrated luminosity J f(E) dE to the photospheric lu- 
minosity from I A. 6l 



A.2. Fitting formulas 

The dependency of the photospheric bolometric luminosity on 
(To is well modeled by a SBPL model (see fig. 3) of the form 

where 

^ = 61n(cro/cTo,c,) (A.7) 

and (To.cr = 42(L52(eQ)3) ^ (defined in eq. dS)). The best fit 
for Al is 

Al « 7 • lO^^\sD.)'^^-^Ll-^ erg ■ s"' . (A. 8) 

The rather sharp break means that for ctq > cro ci the photo- 
spheric luminosity is well represented by (see also eq. ( IA.4b ) 

Lph = 7 • 10'°(eQ)'f L^2Vo,'2^ erg-s-^ ((t xtq) (A.9) 



